Photonic Band Gaps in 3D Network Structures with Short-range Order 
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We present a systematic study of photonic band gaps (PBGs) in three-dimensional (3D) photonic 
amorphous structures (PAS) with short-range order. From calculations of the density of optical 
states (DOS) for PAS with different topologies, we find that tetrahedrally connected dielectric 
networks produce the largest isotropic PBGs. Local uniformity and tetrahedral order are essential to 
the formation of PBGs in PAS, in addition to short-range geometric order. This work demonstrates 
that it is possible to create broad, isotropic PBGs for vector light fields in 3D PAS without long-range 
order. 

PACS numbers: 42.70. Qs, 42.25.Fx, 61.43.-j, 78.20.Bh 



I. INTRODUCTION 

A photonic band gap (PBG) describes a frequency 
range within which light propagation is prohibited due to 
depletion of optical states. The most known structures 
having PBGs are photonic crystals (PhCs) with periodic 
modulations of dielectric constant Since PhCs are 
anisotropic, PBGs vary with directions. To have a 
complete PBG, the gaps in all directions must overlap 
in frequency. This condition is difficult to achieve for 
many PhCs, for example, simple cubic lattices. It is 
therefore easier to produce complete PBGs in more 
isotropic structures, e.g. photonic quasicrystals that 
possess higher rotational symmetry (but no translational 
symmetry) @,||. Photonic amorphous structures (PAS) 
are most isotropic, due to the absence of long-range 
translational or rotational order. Recent studies demon- 
strate that PBGs can be formed in two-dimensional 
(2D) and three-dimensional (3D) PAS with short-range 
order P4Io|. However, the exact physical mechanism 
or condition for the PBG formation in PAS is not well 
understood. An improved fundamental understanding 
of PBG formation would allow researchers to design 
photonic amorphous materials with optimized and 
tunable PBGs. 

In addition to geometric order, structural topology 
plays an important role in forming a PBG. For the com- 
posite dielectric materials consisting of two components 
with different refractive indices, there are two cases 
regarding the topology of the high-index component, 
(i) Cermet topology: the high-index material consists 
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of isolated inclusions, each of which is completely 
surrounded by the low-index material, (ii) Network 
topology: the high-index material is connected and 
forms a continuous network running through the whole 
composite. Previous studies of periodic structures have 
indicated that the cermet topology is more favorable for 
the PBG formation of a scalar wave, while the network 
topology for a vector field Such conclusions also 

apply to PAS. For example, in 2D PAS, PBGs for the 
transverse magnetic (TM) polarization (electric field out 
of plane) are easily obtained with isolated islands of 
high-index materials, because the electric field has same 
polarization direction everywhere and can be regarded 
as a scalar wave. For the transverse electric (TE) 
polarization (electric field in plane) , the electric field has 
varying polarization direction and behaves like a vector 
field, thus it is easier to produce PBGs in connected 
dielectric networks [l2T |. It has been proposed that a 
hybrid structure with a mixture of both topologies can 
possess a full PBG for both TE and TM polarizations 0. 



It is much more difficult to form complete PBGs in 
3D structures. Substantial reductions in the density 
of optical states (DOS) have been demonstrated in 
PAS composed of randomly packed dielectric spheres of 
uniform size 0], as a result of evanescent coupling of 
the Mie resonances of individual spheres. Dielectric net- 
work structures, for example, the photonic amorphous 
diamond (PAD), exhibit much stronger depletion of the 
DOS @, Q. It was conjectured that the tetrahedral 
bonding configuration in the PAD plays an important 
role in the formation of isotropic PBG. However, 
the PAD is constructed from a "continuous-random- 
network" (CRN) ori gina lly developed for modeling of 
amorphous Si or Ge [ljj, thus it is difficult to separate 
the relative contributions of tetrahedral bonding and 
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local geometric order to the PBG formation. Identifying 
the key parameters that determine when a PBG will 
form in PAS is important not only for developing novel 
photonic glasses [lj], but also for understanding color 
generation in nature [l^]. Both cermet and network 
topologies have been found in color-producing PAS of 
many animal species (l6l [l7j . It is also conjectured 
that pseudo PBGs may be formed and responsible for 
non-iridescent coloration of many PAS [l8| . 

In this article, we present a detailed numerical study 
of the DOS and PBGs in 3D PAS. We vary the topol- 
ogy, short-range geometric order, refractive index con- 
trast, and filling fraction to maximize the depletion of 
DOS and the strength of PBG in the absence of long- 
range structural order. This study allows us to identify 
the essential elements for the formation of PBGs in PAS. 
We organize the manuscript as follows. In Sec. [TT1 we de- 
scribe the methods used to generate PAS numerically and 
analyze their structural and topological properties. The 
calculated DOS for both cermet and network topologies 
are presented in Sec. IIII1 together with interpretations 
of the results. In Sec. IIV1 we explore dielectric networks 
with different degrees of structural order to maximize the 
reduction in the DOS. We then conclude in Sec. E] 



II. STRUCTURE GENERATION AND 
CHARACTERIZATION 

A. Spheres packings 

We first study dielectric composites with the cer- 
met topology — high-index dielectric spheres embedded 
within a low-index host material (air). We employ a 
two-stage numerical protocol to generate 'just-touching', 
jammed sphere packings in a cubic simulation cell with 
varying positional order [l9|, First, liquid states 

of monodisperse spheres are cooled at fixed packing 
fraction 4> — 0.60 from an initial high temperature To 
to zero temperature at different rates. In the second 
step, each zero-temperature configuration is compressed 
in steps of A(j> = 10~ 3 followed by minimization of the 
total energy until a static packing with infinitesimal 
particle overlaps is obtained. By varying the cooling 
rate, we are able to create static packings with a range 
of positional order and packing fractions from random 
close packing at 4* = 0.64 to the face centered cubic 
(FCC) structure at <fi = 0.74. In general, the slowly 
cooled samples can be compressed to higher packing 
fractions. Figure [T] (a) shows a cluster of 50 spheres 
from the interior of a jammed sphere packing containing 
1000 spheres at <fr = 0.64. For comparison, we generate 
completely disordered configurations by placing spheres 
randomly in the cubic box with no overlaps at <p = 0.35. 



B. Dielectric network 

We also generate structures with network topologies, 
where the high-index dielectric material forms the 
continuous network, using two methods. For the first 
method, we invert the cermet structure of jammed 
dielectric spheres in air. The inverse structure consists 
of low-index (air) spherical inclusions in a continuous 
high-index dielectric network. By adjusting the radius 
R of the spheres (but fixing their positions), we can 
vary the air fraction 7 in the inverse structure. An 
inverse structure with 7 = 0.8 is shown in Fig. [T] 
(b). At this 7, adjacent air spheres begin to overlap 
and the dielectric material exhibits an irregular topology. 

The second method, which is based on an algorithm 
described in Ref. [H, [2l|, produces more uniform net- 
work topologies than those from the first method. In 
this method, a 3D Delaunay tessellation is performed 
on the sphere centers from the cermet structures in 
Sec. Ill Al Each tetrahedron of the tessellation has four 
facets shared with four neighbors. We then calculate 
the center of mass of each tetrahedron, and connect the 
centers of mass of nearest neighbors by a dielectric rod. 
This creates a tetrahedrally connected dielectric network, 
where each junction (vertex) has four dielectric bonds. 
All dielectric rods have same radius W, but different 
lengths d. By changing W, we can vary the air frac- 
tion 7. A tetrahedral network with 7 = 0.8 is shown in 
Figffl(c). 




FIG. 1. Three examples of photonic amorphous structures: 
(a) jammed packing of dielectric spheres at (f> = 0.64, (b) 
inverse structure of (a) with air fraction 7 = 0.8, and (c) 
tetrahedral network of dielectric rods with 7 = 0.8 obtained 
from the Delaunay tessellation of (a). 



C. Structural characterization 

We now calculate the density autocorrelation function 
and spatial Fourier spectra of the cermet and network 
structures described above. Since the dielectric spheres 
embedded in air and the corresponding inverted structure 
possess identical geometrical properties, we focus only on 
the air spheres and tetrahedral network structures below. 

As shown in the inset to Fig. [5] (a), the 3D spatial 
Fourier transform of the tetrahedral network structures 



3 



displays concentric spherical shells without discrete 
Bragg peaks, which reflects structural isotropy and a 
lack of long-range order. The radii of the shells provides 
the characteristic spatial modulation frequencies of 
the structures. Similar results are obtained for the 
tetrahedral networks generated from the jammed sphere 
packings. The angle- averaged power spectra for both 
sphere and network structures are plotted in Fig. [5] 

(a) . The main peak represents the dominant spatial 
frequency, and its width is inversely proportional to the 
average size of ordered domains (22|. The sphere and 
network structures have similar peak widths, and thus 
comparable domain sizes. 

We also calculated the real-space density autocorre- 
lation function C(Ar) averaged over all angles for the 
sphere and network structures [22| . As shown in Fig. [2] 

(b) , both structures display highly damped oscillations 
of C(Ar). The first peak away from Ar — is located 
at the average spacing a between nearest neighbors. We 
find that the amplitudes of the oscillatory peaks decay 
exponentially [inset to Fig. [2] (b)] with a decay length 
(excluding the first peak) £ r » 0.9a for the sphere pack- 
ings and 1.1a for the tetrahedral networks. Hence, there 
are weak spatial correlations and short-range order in 
these PAS. 
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FIG. 2. Structural characterization of photonic amorphous 
structures, (a) Angle- averaged power spectra of the spa- 
tially Fourier transformed density for jammed sphere pack- 
ings (dashed line) and tetrahedral networks (solid line) versus 
qa/2-K, where q is the spatial frequency and a is the mean spac- 
ing between spheres. The inset shows a cross-section of the 
3D power spectrum for the tetrahedral network, (b) Angle- 
averaged density autocorrelation for the sphere packing and 
network structures. The inset shows the amplitudes of the 
oscillatory peaks of C(Ar) for sphere packings (circles) and 
tetrahedral networks (crosses). 



III. DOS OF PAS WITH CERMET AND 
NETWORK TOPOLOGIES 

In this section, we describe calculations of the DOS for 
jammed dielectric spheres in air, the inverse structure, 
and the tetrahedral networks of dielectric rods using the 
order- TV method [23| . We choose a cubic supercell with 
size 8.7a containing 1000 spheres and refractive indices 



n = 3.6 and 1 for the high- and low-index materials, re- 
spectively. We find that the optimal air fraction that 
yields the largest reduction of the DOS is 7 = 0.75 for 
the dielectric sphere packings and 0.80 for both the in- 
verse structure and tetrahedral network. The DOS was 
ensemble-averaged over five distinct configurations at the 
optimal 7 for each topology, and then normalized by 
the DOS of a "homogeneous" medium with the same 7. 
The latter structure is generated by placing cubic dielec- 
tric voxels (with lateral dimension 0.043a, which is much 
smaller than the wavelength of light A) randomly in the 
supercell. 
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FIG. 3. DOS for (a) jammed dielectric spheres in air with 7 = 
0.75, (b) inverted structures with 7 = 0.8, and (c) tetrahedral 
networks with 7 = 0.8. The wavelength A is normalized by 
the mean spacing between spheres a (average bond length d) 
on the top (bottom) scale. 

As shown in Fig. [31 the maximal DOS reduction occurs 
in the tetrahedral network structure, which is two orders 
of magnitude larger than that for the dielectric spheres 
and inverse structures. For the tetrahedral networks, 
the PBG is formed at normalized frequency d/X « 0.22, 
where d is the average length of dielectric rods and 
d/a = 0.39. The width of the PBG normalized by the 
gap center frequency is ~ 5.5%. The modest reduction 
in the DOS at a/A w 0.41 for the dielectric spheres 
stems from Mie resonances of individual spheres [7| . The 
uniformity of the dielectric spheres allows the coupling 
of their Mie resonances, which of the lowest order for 
isolated dielectric spheres in air occurs at a/A w 0.41. 
In contrast, the air sphere structures have only a small 
reduction of the DOS in the frequency range where the 
tetrahedral networks show a pronounced PBG, despite 
the fact that both structures have dielectric network 
topology and similar degree of spatial correlation. It is 
clear that the dramatic difference in the DOS cannot be 
explained by the small differences in spatial correlations. 

Our studies of jammed dielectric sphere packings 
show that uniformity in the size of dielectric spheres 
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leads to strong coupling of Mie resonances that result 
in a depletion of the DOS. In the inverse structure of 
air spheres, the basic scattering unit is the dielectric 
filling between air spheres. For the tetrahedral network 
structure, the basic scattering unit is centered at each 
junction where four dielectric rods meet. Note that 
in the network topology, the adjacent scattering units 
are connected, in contrast with the cermet topology. 
To compare the uniformity of local scattering units in 
dielectric networks, we calculate the average refractive 
index near the center of each unit. For the tetrahedral 
network structure, we calculate the mean refractive in- 
dex n within a sphere of radius r whose center coincides 
with the center of each junction. We then compute the 
average (n(r)) and its variance V{r) over all junctions. 
For the air spheres, the dielectric junction center is set 
at the center of refractive index distribution within each 
tetrahedron obtained from the 3D Delaunay tessellation 
of the sphere centers. Similarly, we calculate the mean 
refractive index n around each junction center, (n(r)}, 
and V(r) averaged over all junctions. 

In Fig. [4] (a) we show that on average the tetrahedral 
network and air spheres structure have similar distribu- 
tions of the mean refractive index (n(r)) around each 
dielectric junction. In addition, the average refractive in- 
dex for both networks approaches the same value at large 
r since the air fraction 7 is the same for both structures. 
However, the variance V(r) of n for the two network 
structures shows marked differences for all r as shown 
in Fig. 0Jb). The tetrahedral network possesses much 
smaller fluctuations in n from one junction to another. 
Thus, the scattering units are much more uniform for the 
tetrahedral network than those in the air spheres. The 
uniformity of local refractive index distribution ensures 
similar scattering characteristic of individual dielectric 
junctions and facilitates their coupling which leads to a 
dramatic depletion of the DOS. 



tractive index of the dielectric material n. In Fig. [5j we 
show the variation of the PBG for two values of 7 and n. 
Reducing the air fraction below 0.8 leads to a decrease in 
the PBG. A reduction in 7 increases the average refrac- 
tive index of the structure, thus reducing the ratio of the 
index difference (n — 1) to the average refractive index. 
It leads to a decrease of the overall scattering strength. 
In contrast, if 7 is increased above 0.8, there is an in- 
sufficient amount of high-index material to scatter light. 
Thus, an optimal 7 that depends on the refractive index 
n exists at which the scattering strength is maximal and 
the PBG is the largest. As shown in Fig. 0(c), when the 
refractive index of the dielectric material is reduced from 
n = 3.6 to n = 2.8, the maximal DOS reduction shifts to 
7 = 0.72. Moreover, the reduction of the DOS at n = 2.8 
is nearly two orders of magnitude smaller than that at 
n = 3.6. Thus, a large refractive index contrast is needed 
for strong PBGs in photonic amorphous structures. 
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FIG. 5. DOS of tetrahedral networks for two values of the air 
fraction 7 and refractive index n. (a) n = 3.6 and 7 = 0.72, 
(b) n = 3.6 and 7 = 0.6, and (c) n = 2.8 and 7 = 0.72. 




FIG. 4. Uniformity of the local scattering environment for 
the dielectric networks of tetrahedral bonding (solid line) and 
air spheres (dashed line), (a) Mean index of refraction (n(r)) 
and (b) variance V(r) within a distance r from the dielec- 
tric junction center, r is normalized by the mean spacing of 
spheres a (average bond length d) on the top (bottom) scales. 

The formation of a PBG in the tetrahedral network 
structure also depends on the air fraction 7 and the re- 



IV. EFFECT OF SHORT-RANGE ORDER 

In addition to the factors studied above, short-range 
positional order and tetrahedral bond order play impor- 
tant roles in the formation of PBGs in PAS. In this sec- 
tion, we focus on the dielectric network of tetrahedral 
bonding, which yields the largest PBGs, and vary the 
amount of positional and tetrahedral bond order. In par- 
ticular, we tune the positional order of the original sphere 
packings from which the tetrahedral networks are formed. 
The degree of positional order increases with the volume 
fraction of spheres <fi, which varies from 0.35 to 0.69. We 
label the tetrahedral networks (Fig. [B^a)-(c)) generated 
from the sphere packings at <f> = 0.35, 0.64, and 0.69 as 
A, B, and C. 2D cross-sections of the 3D spatial Fourier 
spectra for these structures are presented in Fig. |B] (d)- 
(f). The power spectra of networks A and B consist of 
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concentric shells, but the shell width is notably larger for 
A. Thus both A and B are isotropic structures, but B 
possesses more positional order than A. In contrast to 
A and B, network C features discrete diffraction peaks 
in the Fourier spectrum, and the structure is no longer 
isotropic. 




FIG. 6. Tetrahedral dielectric networks generated from sphere 
packings with packing fraction (a) <p = 0.35, (b) 0.64, and (c) 
0.69. 2D cross-sections of the 3D spatial Fourier spectra of 
the corresponding tetrahedral networks are shown in (d), (e), 
and (f). 

In Fig. [7J we compare the DOS of the tetrahedral net- 
works A, B, and C, with the refractive index of the di- 
electric rods set to n = 3.6. By adjusting the dielectric 
rod radius W, we find that the optimal air fraction for 
all three structures is 7 = 0.8. As expected, network 
A, with the least positional order, possesses the small- 
est depletion in the DOS. However, network C with the 
strongest degree of positional order has a smaller DOS 
depletion than network B. This result contrasts with re- 
cent findings for 2D PAS with air cylinders embedded 
in dielectric materials that show increasing positional or- 
der leads to stronger DOS depletion (22|. To understand 
these results, we must also compare the uniformity of 
the local refractive index distribution and the structural 
topology of the three network structures at fixed radius 
W of the dielectric rods. We find that networks B and C 
have comparable fluctuations in n over all the junctions. 
Thus, local uniformity does not explain the difference in 
the depletion of the DOS for networks B and C. 

To investigate the effects of local topology on the de- 
pletion of the DOS, we compute the tetrahedral order 
parameter @, [9[ 
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where ipjk is the angle between two dielectric rods joined 
at a junction in the tetrahedral network [24]. For a pe- 
riodic diamond network, ipjk — 109.5°, cos(ipjk) = —1/3 
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FIG. 7. The DOS for three tetrahedral dielectric networks (a) 
A, (b) B, and (c) C with positional order increasing from A 
to C. 



for all j and k, and thus £ = 1 at each junction. If the di- 
electric rods are randomly orientated, (Q — 0. In Fig. 
we plot the distributions of ipjk and C for the A, B, and 
C networks, and provide the mean values (ipjk or £), and 
standard deviations and sq. 
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FIG. 8. Characterization of the local topology for networks 
A, B and C. (a) The distribution of angles ipjk between 
dielectric rods j and k at each tetrahedral junction. The ver- 
tical dashed line indicates the angle for the periodic diamond 
structure, ipjk = 109.5°. (b) Distribution of the tetrahedral 
order parameters £ at each junction. The average ipjk and £ 
and standard deviations and s c are also provided. 
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Network A possesses the widest distributions for 
both ipjk and £, which indicates that the local topology 
varies significantly from one junction to another and 
the bond angles within each junction are not uniform. 
The distributions of tpjk and £ are narrower for network 
B, and are peaked at ipjk = 114° and ( — 0.95, which 
indicates that most of the junctions have a similar 
topology to that in a diamond lattice. In contrast, 
network C displays multi-modal distributions for ipjk 
and For example, the £ distribution possesses peaks at 
C = 0.95, 0.72, and 0.5. The first peak reveals that there 
are many junctions with strong tetrahedral order, while 
the second and third peaks reflect the existence of many 
"defect" junctions with low £. Such defect junctions 
are likely located at domain boundaries, and introduce 
irregularity in the local configuration of scattering units. 
Figures [7] and [8] show that photonic amorphous networks 
with strong tetrahedral order and few defect junctions 
have broad PBGs. 



V. CONCLUSION 

In this paper, we calculate the DOS in 3D pho- 
tonic amorphous structures with cermet and network 
topologies. We find that interconnected networks of 
high-index material with uniform dielectric junctions 



and tetrahedral bonding give rise to large isotropic 
PBGs. Further, reduced fluctuations in the refractive 
index around each junction and strong tetrahedral order 
for the angles between the dielectric rods that form 
the junctions enhance isotropic PBGs. High refractive 
index contrast and a low fraction of high-index material 
are also important to PBG formation. We have thus 
identified several parameters that can be tuned to create 
broad isotropic PBGs in photonic amorphous structures 
in the absence of long-range structural order. 
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